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Abstract 

The proposed earlier relativistic mean-field model with hadron masses and coupling 
constants depending on the cj-meson field is generalized to finite temperatures. 
Within this approach we simulate the in-medium behavior of the hadron masses 
motivated by the Brown-Rho scaling. The high-lying baryon resonances and boson 
excitations as well as excitations of the a, uj and p fields interacting via mean fields 
are incorporated into this scheme. Thermodynamic properties of hot and dense 
hadronic matter are elaborated with the constructed equation of state. Even at 
zero baryon density, effective masses of a-uj-p-N excitations abruptly drop down 
for T ^ 170 MeV and reach zero at a critical temperature T 210 MeV. 

Below Tea (at T ~ 190 MeV) the specific heat gets a peak like at crossover. We 
demonstrate that our EoS can be matched with that computed on the lattice for high 
temperatures provided the baryon resonance couplings with nucleon are partially 
suppressed. In this case the quark liquid would masquerade as the hadron one. The 
model is applied to the description of heavy ion collisions in a broad collision energy 
range. It might be especially helpful for studying phase diagram in the region near 
possible phase transitions. 



1 Introduction 



The investigation of thermodynamic properties and phase structure of strongly 
interacting nuclear matter at high baryon densities and temperatures has re- 
cently become important in view of plans to construct new accelerator facil- 
ities (FAIR) at GSI Darmstadt for covering the (5-35) AGeV heavy-ion qyi- 
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ergy range [T]. Increasing interest to this energy region is emphasized by the 
recent proposal for a low-energy campaign at RHIC aimed at identification 
of the critical end-point [2] and by current discussions about the feasibility 
for searching the quark-hadron mixed phase at the Nuclotron-based collider 
(JINR, Dubna) [3]. 

Theoretical predictions for critical baryon densities and temperatures of phase 
transitions depend sensitively on the equation of state (EoS) of both the 
hadronic matter and the quark-gluon matter at high densities and temper- 
atures in the nonperturbative regime. Here we focus on studying the hadronic 
EoS. An EoS of hadronic matter should satisfy experimental information ex- 
tracted from the description of global characteristics of atomic nuclei such as 
the saturation density, binding energy per particle, compressibility, asymme- 
try energy and some other. Some constraints on hadronic models of EoS follow 
from analysis of elliptical flow and data in heavy ion collisions (HIC) . In 
addition to these constraints astrophysical bounds on the high-density behav- 
ior of /3-equilibrium neutron star matter were applied in the recent paper [1]. 

Obviously relativistic effects are important under extreme conditions of high 
densities and temperatures. Microscopically based approaches, as the Dirac- 
Brueckner-Hartree-Fock method, see [5], are very promising but need rather 
involved calculations. Uncertainties are getting larger with increase of the 
baryon density and temperature. E.g., three-body forces are to large extend 
unconstrained. Their poorly-known isospin and temperature dependencies in- 
troduce additional sources of uncertainty. In the quark-meson-coupling models 
the nuclear system is represented as a collection of quark bags. The interac- 
tions are generated by the exchange of a-, uj- and p- mesons treated on the 
mean-field level with quarks, for recent review see [6] . These models might be 
especially useful for the study of the quark deconfinement phase transition. 
But they involve uncertainties due to phenomenological description of the 
nucleon structure. As a more economical approach for describing global prop- 
erties of nuclear matter at the above extreme conditions the relativistic mean 
field (RMF) approach is often used where baryons interact with a, uj and p 
mean fields. Parameters of the model are extracted from the comparison with 
the experimental data. First RMF models used a minimal coupling of nucleons 
with (T, UJ and p mesons [7]. However it proved to be insufficient to appropri- 
ately describe experimental data. Therefore non-linear self-interactions of the 
cr meson have been introduced, see [8]. This approach was latter extended 
to meson fields [9]. As an alternative, RMF models with density dependent 
nucleon-meson couplings were developed (TU]. They allow a more fiexible de- 
scription of the medium dependence. Many other extensions of the models 
have been considered, e.g. models including SU{?>) symmetry [11]. It is not 
our goal here to study and compare different models with each other. Most of 
the models were developed in order to describe a specific domain of nuclear 
physics. Their validity in other regions of nuclear physics either was not con- 
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sidered or they failed to describe them. In Ref. [1] a general testing scheme was 
developed to apply the models to all known nuclear systems. It has been shown 
that among other phenomenological models studied there, the RMF model of 
the EoS suggested in [12] proved to be one of the most efficient model satisfy- 
ing majority of the existing constraints at zero temperature (see KVOR EoS 
in Table 5 in Ref. [1]). Therefore focusing on the further application to HIC, 
the given model will be generalized here to finite temperatures. 

Following Ref. [T2| we assume a relevance of the (partial) chiral symmetry 
restoration at high baryon densities and/or temperatures [13] manifesting in 
form of the Brown- Rho scaling hypothesis [H]: Masses and coupling constants 
of all hadrons decrease with a density increase in approximately the same way. 
Note that most of the models use the constant a, u;, p effective masses. Some 
models introduce field interaction terms leading to an increase of the a, cu, 
p effective masses with the increasing nucleon density, e.g. [TnHS]. Contrary, 
Ref. [12] follows the Brown-Rho scaling hypothesis and scales the quadratic 
(mass) terms of a, uj and p fields as well as the baryon mass by a universal 
scaling function $. The scaling function $ is assumed to be dependent on the a 
mean-field. This provides thermodynamical consistency of the model. In order 
to obtain a reasonable EoS, the meson-nucleon coupling constants should be 
also scaled with the a mean field. Differences in the scaling functions for the 
effective masses of a;- and p-fields and their couplings to a nucleon allow one to 
get an appropriate density-dependent behavior of both the total energy and 
the nuclear asymmetry energy, in agreement with constrains obtained from 
measurements of neutron star masses and surface temperatures [1]. 

Our main goal here is to construct some effective model of EoS that would 
incorporate the decrease of hadron masses and couplings with increase of the 
baryon density ub and temperature T and, simultaneously, would fulfill vari- 
ous constraints known from analysis of atomic nuclei, neutron stars and HIC. 
Our consideration is based on the generalization of the KVOR model [1] to 
finite temperatures. We aim to test its suitability to description of properties 
of hot and dense matter formed in HIC. Besides the nucleon and meson mean 
fields we include low- lying non-strange and strange baryon resonances (A^, A, 
A, S, S, S*, S*, and fi), meson excitations o"(600), p(770), 0^(782) constructed 
on the ground of mean fields , and the (quasi) Goldstone excitations 7r(138), 
-R'(495), 77(547). We add here their high mass partners in the SU(3) multi- 
plet i^*(892), rj (958) and (/?(1020). All corresponding antiparticles are also 
included. All states are treated within quasiparticle approximation. 

We restrict ourselves by taking into account only the large N^. ground states 
of baryons and mesons [T7]. Generalization to higher mass resonances is, of 
course, straightforward. However we will drop them from our consideration 
being guided by the following arguments: The next in mass state that does 
not enter the multiplet is the A(1405) state. It manifests in the kaon scattering 
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data as a quasimolecular state near the kaon-nucleon threshold. However in 
description of a broad kaon energy region under discussion in this work the 
A(1405) hyperon does not manifest as a pole-hke term [IB]. Ref. [T^ presented 
arguments that A(1405) should dissolve in matter. Thereby we will exclude 
A(1405) from our study. Moreover as it is conjectured in Ref. [20], higher mass 
baryon resonances can be understood as composite particles. Their widths in 
matter are expected to be quite large. Ref. [20] reproduces particle scattering 
data assuming that the lowest baryon octet and decouplet are only relevant 
degrees of freedom. Therefore we do not incorporate the higher resonances 
within our quasiparticle model. Besides, the higher mass particles are con- 
sidered, the less they contribute to thermodynamics, and the less one knows 
about their interactions. Following above argumentation and in order not to 
complicate consideration by introducing dependencies on unknown parameters 
we accordingly cut the baryon particle set. 

Although free a-, p-, tu-mesons are rather heavy, their effective masses in 
matter essentially decrease. Therefore we include excitations of these fields as 
well. Sometimes the (5[ao(980)]-meson is incorporated in the RMF scheme |21j . 
Its role in RMF models is similar to that of the p-meson except that inclusion 
of 5[ao(980)]-meson allows for mass splitting between proton and neutron. 
Since 5[ao(980)] coupling constants to other particles are unknown and the 
mass is larger than that for uo we do not incorporate 5 in our scheme. 

To construct a practical description, the particle interaction with a-, uj- and 
p-meson fields is treated only in the mean-field approximation. The fermion- 
fermion hole loop diagrams for boson propagators and the boson-fermion loop 
diagrams for fermion propagators are disregarded. Thereby we omit the p- 
wave pion-baryon and kaon-baryon interaction effects, though these effects 
are important in high-baryon density regime [T5]|22|23] . At high temperatures 
the fermion-antifermion loops in boson propagation and fermion-boson loops 
in boson propagation might become very important [21] . We however postpone 
study of these effects to a future work. At sufficiently large baryon densities 
there might appear condensates of some (quasi) Goldstone bosons, like K~ and 
K^, p~, and might be rj. To deal with stable ground state we include a self- 
interaction of (quasi) Goldstone boson fields. With this model we construct the 
EoS as function of the temperature and the baryon density and apply it in a 
broad density-temperature region. Below for brevity we call thus constructed 
model as the scaled hadron mass and couplings (SHMC) model. 

The paper is organized as follows. In sect. [2] we formulate the Lagrangian of 
the model. Then in sect. |3] the energy density for a system at finite baryon 
density and temperature is constructed. In sect. |4] parameters of the model 
are determined by fitting them to available data at zero temperature or by 
exploiting the symmetry relations in case when there are no data. Sects. [5] and 
E] are devoted to evaluation of thermodynamic properties of the constructed 
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EoS for T = and T 7^ 0, respectively. In sect. [7] we apply our model to 
HIC. Some concluding remarks and perspectives are given in sect. [H Lengthy 
formulae for different terms of the energy density and a scheme how to cal- 
culate condensates, when they occur in some (T, n^) range, are deferred to 
Appendices A and B. Throughout the paper we use units h = c = 1. 



2 Lagrangian 



Within our model we present the Lagrangian density of the hadronic matter 
as the sum of several terms: 



Let US describe each term in ([1]). 

The Lagrangian density for baryons interacting via a, u, p is as follows, cf. [12 



>Cbar = E [^b{tD- 7) - ml ^fe + t^QH^V] . (2) 

be{b} 

The long derivative D is given as follows 



D^ = d^- iV t^Q 6^0 + i g^b Xu^ + iQpb Xp P/. , (3) 

where g^b and gpb are coupling constants and Xuj, Xp are coupling scaling 
functions which will be determined below, = {ujq,uj) and = (Po'P") are 
u- and p-fields with a = 1, 2, 3. Here 7^ are Dirac matrices, Ub is the particle 
(or antiparticle) density of baryon species b,V = — e$ei is the electric potential 
$01 measured in the electron charge units. 

The baryon set {b} that we use is presented in Table 1. In reality masses of 
charged and neutral particles of the given species are slightly different. We 
ignore this difference that allows us to use isospin invariance for ^ = 0. 

The a-, u-, p-meson contributions to the Lagrangian density 



vCmes = E '^rn, with {m} = a,uj,p (4) 
render, respectively: 
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Table 1 



Properties of the basic baryon set (masses are given in MeV) 
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1193 


1193 


1193 
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1 





-1 
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-1 





1 





-1 
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-1 


-1 


-1 


Sb 


1 
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2 


1 
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1 
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1 

2 



continuation 
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1 
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^b 





-1 
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-1 





-1 


-1 


^b 


-2 


-2 


-1 


-1 


-1 


-2 


-2 


-3 


Sb 
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2 
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2 


2 



H ^ , ^^/i^ 



(5) 



^P'' = 2 2 ' ^ 

V = [{td' -V + g,Xp ) [(^5° -V + g,Xp P^o^) Pck]^ 

- iVp^ftP -m;2|p- |2, p-^ = (pi -ip2)/V2. (7) 

Here Qp coupling is responsible for the self-interaction of the charged and 
neutral species, Xp is the scaling function. Non-Abelian p-p interaction with 
dp = QpN is motivated by the hidden local symmetry approach, cf. [25], where 
p-meson is introduced as a non-Abelian gauge boson. Nevertheless this pos- 
sibility is often disregarded and one uses the simplest form with Qp = cf. 
[26] . In a sufficiently dense asymmetric nuclear matter the presence of the 
self-interaction may result in appearance of the charged p-meson condensate 
characterized by non-zero p~ mean field instead of pO one, cf. [T21I27] . 

Following [12] we use the cr-field dependent effective masses of baryons 
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ml/mb = ^biXaCr) = I - Qab Xa cr / nib , he {h} (8) 

with the baryon set {6} defined in Table 1 and mass terms of the mean fields 
are 

'^mMm= |$m(XaCr)| , {m} = (T, CJ, p , (9) 

where g„b are afe-coupling constants. We have introduced the absolute value 
of $m(Xo-C") to indicate that only {rn*^Y mass terms, as they appear in the 
Lagrangian, have physical meaning. This observation is important to interpret 
the situation when $m(Xcr cr) becomes negative, see sect. E] below. 

For the sake of simplicity we scale all couplings (7^6 by a single scaling function 
Xo-(c"), and all g^^^i Qpb by Xuj{(^) and Xp{^) scaling functions, respectively. Thus 
all scaling functions depend only on a. The idea behind that is as follows. The 
cr field can be expressed in terms of the ud quark condensate. The change of 
effective hadron masses and couplings is associated namely with modification 
of the quark condensate in matter. Thus we consider the a field as an order 
parameter. The a excitations are then treated as fiuctuations around the mean 
value of the order parameter. Similarly long-scale fiuctuations are treated in 
the Landau phenomenological theory of phase transitions. 

The dimensionless scaling functions $6 and as well as the coupling scaling 
functions Xm depend on the scalar field in the combination Xcr(o") a. Therefore 
for further convenience we introduce the variable 

f = QaN Xa (j/rriN . (10) 

Following [12] we assume an approximate validity of the Brown-Rho scaling 
ansatz in the simplest form 

$ = $^ = $^ = $^ = $^ = 1 _ /. (11) 

We keep the standard form for the non-linear self-interaction term (potential 
U) of RMF models, but now in terms of the new variable /, and using ffTOl) it 
can be rewritten as follows: 

u —mj^y^- J -I- - / J — ^ I ^ . l,i/J 

Two additional parameters, h and c, allow us to accommodate realistic values 
of the nuclear compressibility and the effective nucleon mass at the saturation 
density. An extra attention should be paid to the fact that the coefficient c 
must be positive to deal with the stable ground state. 
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The contribution of the electric field Cc\ to the Lagrangian density is: 



= ^_(VVf . (13) 

Coulomb effects are responsible for a deviation of the low-momentum ir^/n^ 
rates from unity in HIC of isospin-symmetric nuclei, and for some other effects. 
It is important to include the Coulomb term for the description of mixed 
phases in dense neutron star matter, cf. [28||29] . 

There are mean-field solutions of the Lagrangian iZbar+'^mcs+'Cci = J2be{b} '^b+ 
J2m(^{m} l^m + To thesc terms we add the Lagrangian density 



HgoXA = E {9} = vr±'°(138); i?°(495); r/(547). (14) 

9&{9} 

These particles are often treated as (quasi) Goldstone ("Gold") bosons within 
the chiral SU(3) symmetrical models. Therefore we may not to scale their 
masses and couplings, as we have done for {m} = a,u!,p. At rather small 
baryon densities there are no mean-field solutions of equations of motion which 
follow from Ccoid- Such solutions may however arise at sufficiently large baryon 
densities signalizing on condensations of these fields. On the other hand, we 
observe, cf . [12] , that for the case of spatially homogeneous system equations 
for mean fields and thus mean- field solutions do not change if we replace a-, 
LU-, p-fields by the scaled fields x^<^, Xc^^ and XpP provided ^t, = ^ 
and x'p = Xp- If we wish to extend this symmetry to the case when Goldstones 
are included, in addition to scaling of masses we should scale couplings, g^g = 
9mg Xm- Bclow we wiU test both possibilities g*^g = g^g and g*^g = g^g Xm, 
and refer to them as versions without and with scaling, respectively. 

The contribution of the pion Lagrangian density £^ into Eq. (fT4l) is given by 

£^ = £^o + £^±, (15) 

= if ^ 1 m^ = m^- g^^a; (16) 

£^± = (ido -V + g*^uJo + ^;,Po)vr- [{ido -V + g*^ujo + g;yo)^l^ 

-|V7r-|2-m:V"l' • (17) 

The kaon Lagrangian density contributes as 

= [{do - iqxV + iglK^Q + z^^^rsp;])^^^] 

X [(.9o - iqxV + iglK^o + ig*pK^-ipl)K] (18) 
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- WK^VK - m*^K^K, m\ = niK - qIk'J, Qk = ^{1 + , 
with K = [K, K^) , Ta is the Pauh matrix. 

We follow [30] and present the Lagrangian density of the t] meson as 

+ ^ E d'^V d,v- (19) 

Here /tt = 93 MeV is the pion decay constant. For nucleons the constants 
S^AT and KriN are estimated from the scattering data: S^tv — 280 ± 130 MeV 
and KriN — 0.40 ± 0.08 fm. Other parameters in (fT9!) are not known. Ref. [30] 
used a large value for the i^'A^-sigma term, T,kn- Below we argue for a smaller 
value of Tjkn- Therefore we will test a sensitivity of the //-description to the 
coupling variations. 

When the total Lagrangian is known, one can derive equations of motion for 
every field. Even for low baryon density, equations of motion for a, uo and p 
and V allow mean-field solutions do, ujq, Pq, and Vq. Therefore we use: 

cr = (Jo; cj^ = cjo p"^ = Ro Sa3 S^o; V = Vq. (20) 

Only for isotopically asymmetric matter (A^ ^ Z) we have i?o 7^ 0. As we 
have mentioned, if the baryon density increases above a critical value and 
N ^ Z, there may appear another solution with i?o = but with a non-zero 
solution for the charged p-meson mean field, p~j^ 7^ 0, cf. [T2|27| . For the sake 
of simplicity we disregard such a possibility in the present work. 

Similarly to the case of self-interacting p meson fields, there exist higher or- 
der terms in Cqo\a describing self-interaction of the fields. Using approximate 
SU(3) theory these terms can be presented as 

^Goid = A E(<?') V4; $ = (vTi, 7r2, TTg; i^i, K^, Ks, K,; 77) (21) 
9 

with a positive self-interaction coupling constant A = const ~ 1 and redefined 
fields 7r± = (vti ± i7r2)/V^, vr^ = vrg, = {Ki ± iK2)/V2, = {Ki + 
iK2)/V2, K° = (Ki - iK2)/V2. Eq. has the simplest form ahhough we 
could use self-interaction terms with different couplings for different particle 
species. 
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The remaining terms SC^cc in the Lagrangian density ([T]) are due to the vector 
mesons K* and (f, and the glueball-hke state rj . One could treat (p like u with 
similar scaling of the effective masses and couplings. Due to a high value of 
strange quark mass one can expect that K* and (p couplings are less than 
those for u. Since little is known about interactions of these particles and not 
to complicate further consideration by introducing unknown parameters we 
consider K* and ip as free particles in present work. To our knowledge there 
is no information about values of rj 6-couplings. Therefore, being conservative, 
we put them zero treating 77 also as a free particle. 

As we have mentioned, condensates of some (quasi) Goldstone fields may ap- 
pear at some specific conditions. In this case their equations of motion acquire 
mean-field solutions, as those for a, uj and p. In such cases for neutral fields 
(the strangeness and electric charge being zero) the stability of the ground 
state is provided only due to presence of the self-interaction, see Eq. fl^T]) . 

To single out quasiparticles (excitations) from mean fields, in the Lagrangian 
£mcs one should do replacements uq = uJq^{<j) + uj', Rq = -Rq (cr) + R^, = 
(f)^\(j) + iu = Lu' and p = p' . Here Uq, Rq, 0'^' are the mean (classical) field 
variables and u'^, (po)^, {p'±Y, <t>' are responsible for new excitations. Then 
we expand the Lagrangian density retaining only quadratic terms in the fields 
of excitations. The coefficients at non- derivative quadratic terms are read as 
squared masses of excitations. We recognize that effective masses of the zero 
and spatial components of vector fields are equal and the gauge conditions 
^^'uJ^, = 0, d'^p''^ = are fulfilled. 

Equations of motion for mean fields and for excitations are obtained by the 
variation of the total action. If mean-field terms are rather large and excitation 
contributions are small, one may disregard the excitation terms in equations of 
motion for the mean fields and neglect self-interactions of excitations. However 
one should keep interactions of excitations with mean fields in equations of 
motion for excitations and in their thermodynamic quantities. 

Since there are no experimental indications of condensation of (quasi) Goldstone 
bosons in the regimes of HIC, we will focus our further discussion on the case 
when condensates do not occur, paying a special attention to situations when 
condensation is possible. 

We checked that minimization of the energy density with respect to the a 
field produces an equation of motion for this field being in agreement with the 
Lagrange equation after its Gibbs averaging. Thus to find the squared effective 
mass of the a excitation (i.e. of the fluctuation of the order parameter), we 
take ci;o(cr), Roicr), plug them in the energy density, put a = a'^^ + a , and find 
the second derivative of the energy density in respect with cr'^\ see Eq. (!66|) in 
Appendix A. 
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3 Energy density at finite temperature 

Let us assume that the system volume is sufficiently large and surface effects 
may be disregarded. Thus only spatially homogeneous RMF solutions of the 
equations of motion are considered. To simplify expressions we will treat all 
quantities in the rest frame. Generalization to the arbitrary moving inertial 
frame is obvious. 

The thermodynamic potential density fl, pressure P, free energy density F, 
energy density E and entropy density are related as 

E = F + TS, F[f, uo, i?o, T] = Y, ^^^n^ + Q = -P, 

i 

_ dF 
drii ' 

Summation index i runs over all particle species; rij are particle densities, see 
Eq. (1^ . Chemical potentials /ij enter Green functions in the standard gauge 
combinations + /ij. 

The energy density can be presented as the sum of the mean a-, u;-, p-field 
contributions as well as contributions of baryons and all meson excitations. So 
we have 

E[f,u,,Ro,T]=Y. E,[f,T]+ E^''[f,uo,Ro,T] 

feg{b} m€{m} 

+ -E'bos.cx. [/, ^o,Ro,T] = Emf + -E'bos.cx. • (24) 

The ffist two sums (resulting in the term -Emf) are included in every RMF 
model but with smaller set {b}, whereas the boson excitation term -Ebos.ex. is 
obtained here beyond the scope of the RMF approximation. 

Although all terms are functions only of / and T, we will present them also as 
functions of lvq and Rq in such a way that the values of the uJo{f) and Ro{f) 
mean fields can be found by minimization of the energy at fixed /. Then 
co'o(/) and Ro{f) are plugged in the energy density functional that becomes 
function of / only. So the equilibrium value of / can be found by subsequent 
minimization of the energy in this field. 

Since Eijos.cx.[f,^o, Ro^T] depends on the mean fields, its minimization pro- 
duces extra terms in the mean-field equations. Within the approximation of 
rarefied gas of excitations used in this work, we treat excitations perturba- 
tively thus omitting these extra terms. Therefore, we assume that -E'bos.part — 
^bos.part[/^'',c^^^,i?^^,T], where f^^,uj^^,R^^ are found by minimization 
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(22) 
(23) 



of the energy E = Emf, i-e. without inclusion of the boson excitation term. 
Thus our equations of motion for mean fields are: 



d 



EMF[f,UJo,T] = , 



d 



EmfI/, Ro,T] = 



(25) 



dR 



-0 



and 



d_ 



EMF[f,u;o{f),Ro{f),T]=0. 



(26) 



At variation of the energy density one should not vary over the particle oc- 
cupation numbers. Note that if the interaction of excitations were included 
within the self-consistent Hartree approximation, in Eqs. ( |25|) . (!26|) . we would 
minimize the total energy E rather than Euf- It would however additionally 
complicate the solution of the problem. We postpone the study of this possible 
model generalization for a future work. 

It is noteworthy that to obtain equations of motion, instead of the energy 
density one could vary the thermodynamic potential Q, cf. [2H]- 

At the resonance peak the vacuum A-isobar mass width is F^'*^ ~ 115 MeV. 
In reality Fa is the temperature-, density- and energy-momentum-dependent 
quantity. For low A-energies the width is much less than F^'^'^. A typical A 
energy is — m*^ ~ T. Thus for low temperature, T < ep {^f '^^ the nu- 
cleon Fermi energy) the effective value of the A-width is significantly less 
than F^'^^. At these temperatures the quasiparticle approximation does not 
work for A's but their contribution to thermodynamic quantities is small. 
When the temperature is > m^r there appears essential temperature contri- 
bution to the width and the resonance becomes broader [31]. A's essentially 
contribute to thermodynamic quantities. Only for temperatures T > T^^{T) 
the quasiparticle approximation becomes a reasonable approximation. 

In reality p- and cx-mesons also have rather broad widths. The observed en- 
hancement of the dilepton production at CERN, in particular in the recent 
NA60 experiment [32] on production, can be explained by significant 

broadening of the p in matter [33] , though decreasing of the p mass could also 
help in explanation of the data |34p1 . Besides, the p width might increase 
with further decrease of its effective mass [32] • Also particles which have no 

^ As demonstrated in [33] the calculated large mass shift is mainly caused by the 
assumed temperature dependence of the in-medium mass. Inclusion of this tem- 
perature dependence modifies the scaling hypothesis originally claimed by Brown 
and Rho. Some arguments on what the proper mass-scaling predicts for dilepton 
production in HIC, e.g. NA60, were given in [35] . 
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widths in vacuum like nucleons acquire the widths in matter due to collisional 
broadening. Their widths grow with the temperature increase, cf. [23]. As we 
have argued above in case with A isobars, the quasiparticle approximation 
may become a reasonable approximation at sufficiently high temperature, if 

T > r(T). 

The problem becomes much more involved, if one tries to treat particle width 
effects consistently. Therefore, to simplify consideration we use the quasipar- 
ticle approximation in the present work for all particle species in the whole 
temperature and baryon density range of our interest. 

Now let us subsequently consider all energy terms in Eq. 
3.1 The baryon contribution 



The contribution of the given baryon species 6 to the energy density is as 
follows 



J 2 

E,[f, T] = {2s, + 1) I ^ (/, + f,)^mf{f)+p^- n, {V + /ieh) 



+ t^ nk{V + fi,^), p=\p\. (27) 

The spin factor s;, = 1/2 for and hyperons, while Sb = 3/2 for the A- 
resonance, see Table 1. The Fermi-particle (baryon) occupations. 



fb = ; , , (28) 

exp[(y^mf + p2_^*)/T] + l 

fb = , , ^ , (29) 

exp[Umf + p^ + fxl)/T] + l 



depend on the gauge-shifted values of the chemical potentials 

Aifo = h /ibar + tl A^str + if^ch + V) - Q^b Xuj^O " ^1 Qpb Xp Rq ■ (30) 

The baryon chemical potential of the b species is fib = tb yWbar; and the corre- 
sponding strangeness term is /i^ = tl figtr- As is seen, the electrical potential 
V ^ V + fich is shifted by the charge chemical potential /Xch related to the 
isospin composition of the system. Suppressing Coulomb effects one drops out 
the shifted value of V. Sometimes instead of fich one introduces the isospin 
chemical potential, cf. [37] . 
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3.2 Mean-field contribution 



It is convenient to introduce the coupling ratios 

Xmb = gmb/gmN, {m} = (J, U, p, (31) 

and, instead of Xm? another variables 

Vm{f) = ^l{f)/Xl{f), (32) 

since the energy density depends namely on such combinations rather than on 
$m and Xm separately. 

Using these new variables the contribution of mean fields to the energy density 
is given as follows: 



Er[f] = '^vAf) + Uif), (33) 



Cl (Ef,g{b}^a;b (^b - rif,))^ 



2C2 



9ujN Xio ^0 



Cl (Efeema^o^b {^h - rib)) 



^nVM) 



(34) 



The net baryon density is defined as follows 



{Ub - Ub) = Ub , 

be{6} 

where the partial baryon and antibaryon densities for the species h are 



(35) 



nb = {2sb + l)J^fb, nb={2sb + l)J^fi 



dp - 



(36) 



Renormalized constants are 



C„ 



fTT-N 9mN 



nrir 



(37) 
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Values of the parameters used will be specified in sect. 4 below. Similarly, for 
the p mean-field contribution we have 



CI 



n 



B) 



m 



N 



Vpif) 



>m 



N 



Vpif) 



2 CI 



QpN Xp Ro 



H 2 



2m% 



ripif)\ 



The isotopic charge density in the baryon sector is given by 



n 



B 



2E 

feG{b} 



rib) Xpb 



(39) 



As one can see, the isovector baryon density plays the role of the source 

(3) 

for the p-meson field Pq = Rq . Therefore for the iso-sjTiimetrical matter 
(A^ = Z) one has nf = and Ef^ = . 

The net density of strange baryons and mesons reads 



"str = Y tl [rib - Ub) - Uk- - n^o - riK*- - n^*o 
be{b} 

+ riK+ + riKO + nK*+ + nK*o 



(40) 



We assume that all strange particles are trapped inside the fireball till the 
freeze-out. Therefore the total strangeness is zero. In this paper we do not con- 
sider the possibility of a mixed phase: Strange clusters surrounded by normal 
matter. This possibility arises since the charge (electric charge, baryon charge, 
strangeness, etc) can be conserved only globally rather than locally |38j . 
We put locally ristr = 0. Then this condition determines the value of the 
strangeness chemical potential p^tv 

Similarly, we may introduce the electric charge density 



riK' 



rich = ^bi^b - rib) + ri^+ + riK+ 
fee{fe} 

- n^- - riK- - Up- - riK*- 

assuming that it is conserved locally. The quantity rich 
the value of the charged chemical potential pch_- 



(41) 

{Z/ A)nB determines 



Our SHMC RMF energy density functional depends on four particular com- 
binations of the functions, rj^^^p^^if) and U{f). Note that the dependence on 
the scaling function 77^ can always be presented as a part of the new potential 



U obtained by means of the replacement U ^ U 



(1 - Vaif)) , and 
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vise versa, so the potential U can be absorbed in the new quantity 7]^. Thus 
actually only three independent functions enter the energy density functional. 
Eq. dMD together with Eqs. (I2ZD> (I331)> (I31D> dM]) demonstrate explicitly equiv- 
alence of mean-field Lagrangians for constant fields with various parameters 
if they correspond to the same functions ?7p,a;(/) and ?7o- (either U{f)), with 
the field / related to the scalar field a through Eq. ffTOl) . 

3.3 Bosomc excitations 

To find the total energy (1241) one should yet define the contribution of bosonic 
excitations. The energy density of boson excitations is the sum of partial 
contributions 

i?bos.cx[/, u;o, Ro, T] = + + E^^^^' + Ef"'* 

+ Ef'' + EP^'' + E^:' + E^r' + E^""'' . (42) 

Explicit expressions for the partial contributions can be found in Appendix A. 

As was mentioned, in the present paper we consider a non-interacting gas 
of excitations. Thus, in order to get E^^^^ we expand EMF[c,^o{'^)y Ro{<^),T] 
in Eq. in 6a = a — a^^. Linear term in 6a does not contribute due to 
subsequent requirement of the energy minimum in a'^^. From the second order 
term we extract 

(mr**)' = d^EuF[a.uj^{a),Ro{a),T]/da\ (43) 

Here particle occupation are not varied as in fl26l) . In the gas approximation 
we may drop the higher order terms in 5a. Effective masses of u and p prove 
to be the same as those follow from the mean-field mass terms 

m^'* = mj<l>^(/)|, mp* = mp|$^(/)| . (44) 

As have been mentioned, the simplifying ansatz (llip is used in present work. 

At certain conditions Bose condensates of some boson species may occur. 
Below we will demonstrate that our choices of couplings do not allow for 
Bose condensates of excitations in the temperature-density region which we 
will discuss in application to HlC However condensates may appear for other 
possible choices of couplings and if a broader density-temperature interval is 
considered. Explanation how to include Bose condensates if they appear is 
given in Appendix B. 
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4 Choice of SHMC model parameters and scaling functions 



Parameters of the RMF model, C^, C^^ Cp, and the self-interaction potential 
U, are to be adjusted to reproduce the nuclear matter properties at the satu- 
ration for T = 0. Usually they are fixed by values of the binding energy ebind, 
nuclear saturation density ub = uq and symmetry energy coefficient asym, 
which are known within some error bars. We will use the same basic input 
parameters as in Ref. jl2] : 

no = 0.16 fm-^ ebmd = -16 MeV, asym(no) = 32 MeV. (45) 



The saturation baryon density and the binding energy are related as 



dE[nB-J] 



dn 



B 



— E[nQ] /(no)] =mN + e^nd 
no 



(46) 



"0 J{"o) 

and the compressibility modulus is given by 



= 9no 





( d^E 


)' 


'd^E 




-1' 






no, f {no) J 


dp 


no,f{no)_ 





(47) 



Here /(no) is a solution of Eq. (l26l) at the density Up = Un = no/2 . The pa- 
rameter Cp is determined from the symmetry energy coefficient of the nuclear 
matter 



"sym 



{ub) 



ub d"^ 
8 dn^ 

8m% r]p 



E {ub - Up, Up) 



+ 



vr^ Ub 



I'M 



np=n/2 



(4J 



Pfn 



Pfn being the nucleon Fermi momentum in iso-symmetrical matter (A^ = Z). 



We use the modified Walecka model with a non-universal scaling of masses and 
couplings (referred as the MW(n.u) model in [12] and as the KVOR model 
in [1]). This model matches the Urbana-Argonne EoS (A18+(5i;+UlX*) [32] 
for the baryon densities below 4no at T = 0, which correctly reproduces the 
maximal neutron star mass M^ax — 2 Mq and gives sufficiently large threshold 
density for the direct Urea reaction n^^^ to be in agreement with the neutron 
star cooling phenomenology [10|. The Urbana-Argone (Al8-|-(5f +UIX*) EoS 
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is derived within a microscopical variational theory of nuclear matter. It em- 
ploys the non-relativistic paired NN potential extracted from the analysis of 
the scattering data, includes the boost u^/c^-order corrections and incorpo- 
rates a three nucleon interaction. However due to using of the non-relativistic 
potential the Urbana-Argone EoS violates causality for ra^ > 4no and T = 0. 
In Ref. gl] the AlS+^t'+UlX* EoS is fitted for n < Auq but the causality 
problem is solved for higher densities. We call this modification as the HHJ 
model. Parameters of our SHMC (KVOR-based) EoS are fitted in such a way 
that the energy (including the symmetry energy) and the pressure are very 
close to those for Al8+5t;+UlX* (and HHJ) for ra^ < 4no and T = both for 
the Z = and N = Z cases. Since our EoS is based on RMF calculations, no 
causality problem arises. 

An appropriate behavior of the EoS for T = is obtained with the scaling 
factors introduced in the same way as in [12] 



Va = — = l, vM) = , ^ f , 49 



with parameter. If one puts z = 0, the standard RMF version without 

a scaling is covered. Effective nucleon mass and the compressibility coefficient 
are 



mlj{no)/mN = 0.805, K = 275 MeV. (51) 

Note that, if we chose a smaller values of m*j^{no)/'mN, we should simulta- 
neously increase K to supply positivity of the c constant in the interaction 
potential U. 

The used here values of other parameters of the SHMC model are the same 
as in the KVOR one 



^ = 0.65 : 
^ = 87.60C 

b = 7.7346 X 10"^ , c = 3.4462 x 10"^ . (52) 



= 87.600 , = 100.64 , Cj = 179.56 

-3 „ Q /I 1 n— 4 



As we will show below, the scaling ratio 77^^ is in the interval 0<1 — ?7(^^0.15 
at T = for all baryon densities of our interest. Then rip is always finite and 
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positive. Considering symmetric nuclear matter N = Z at a. high temperature, 
we will demonstrate that / — > 1 (for T 2^ 210 MeV) that corresponds to 
1 — ?7;^ ~ 0.32. However, already for a smaller value 1 — rj^ the ratio r]p in 
(1501) has the pole. Actually, in the case N = Z the = and the problem 
does not arise. Nevertheless, if we wanted to describe the high-temperature 
regime for N ^ Z in a. similar way as for N = Z,we would need to correct the 
above expression for rjp. Thus we suggest instead of Eq. fl50|l to use its Taylor 
expansion 
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n=0 



2 



(53) 



that solves the pole problem. To reproduce (!50|) in the range < l—t]^ ^ 0.15 
it is sufficient to take 10 terms in (l53l) . Certainly we could introduce other 
parameterizations for the scaling functions. However we will use advantages 
of the KVOR model have been demonstrated in [HIT^ in application to cold 
nucleon matter. Therefore we use the choice fH9|l . (l53l) . 

Like in the KVOR model, we use a rather large value of the Dirac effective 
nucleon mass at the saturation (see Eq. (!5T|) ) as compared to the most of RMF 
models describing finite nuclei. The latter models usually assume m*^/mN 2^ 
0.54 -h 0.7, cf. [12]. Only in this case these models allow one to appropriately 
describe the orbital potential in finite nuclei. Note that if effective meson 
masses are used instead of free ones, the depth of the diffuseness layer of the 
nucleus is changed, affecting the value of the spin-orbit potential. This rises 
hope that the spin-orbit potential problem could be resolved if we solved space- 
inhomogeneous equations in our model. Indeed, already in the framework of 
the standard non-linear RMF model making use of a smaller value of the 
a mass, Ref. [22] permitted to fit the nucleon density profiles. Anyhow the 
description of finite size effects as well as the solution of the mentioned problem 
are beyond the scope of our consideration in the present work. Noteworthy 
that one should distinguish between the Dirac and Landau effective masses. 
The latter value is manifested in the nucleon spectra. Experimental data seem 
to favor the Landau effective nucleon mass being close to the free nucleon 
mass, cf . [^44|45j . Recent calculations of both the Dirac and Landau effective 
masses within the Dirac-Brueckner-Hartree-Fock approach give my/rriN — 
0.7 for the Dirac mass and ~ 0.8 for the Landau mass (see Fig. 2 in [16]) 
whereas pure RMF model calculation and that with the energy-dependent 
correction produce values 0.76 and 0.92, respectively [IS]. Since the ratio 0.92 
for the Landau mass obtained within the model [15] is still low to explain 
nucleon spectra it can be considered as an argument that within this model 
the corresponding ratio my/mM for the Dirac mass should be higher than 
0.76. 
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Now we will demonstrate that the SHMC model describes the nucleon optical 
potential in an optimal way. Optical potentials for a proton or neutron passing 
through the cold (T = 0) nuclear matter are introduced as, cf. [ITPS] : 

U^^, = e-^p^ + ml, N = p,n, (54) 

where e is the nucleon energy and p is the 3-momentum. Substituting from 
equations of motion 

y/p"^ + m*j^ = e - q^n Xco c^o - t% gpN Xp Ro + t^ V , (55) 



U^pt = e - \/(e - gu;N Xu> ^0 - tN9pN Xp Ro + tN^y - - 1), (56) 

t% and are given in Table 1. For the N = Z, V = case, proton and 
neutron optical potentials coincide. 

Energy dependence of the nucleon optical potential for ub = Uq, N = Z 
is shown in Fig. [H Calculated results are presented at T = and there- 
fore coincide with those for the KVOR model. The band is the optical po- 
tential extracted from the data [19] and recalculated to the case of the 
infinite nuclear N = Z matter in [17] . Different lines are RMF calcula- 
tions within the standard Walecka model for various effective nucleon masses 
fn*N/^N = 0.54, 0.71, 0.8, 0.85 [17]. We see that the resulting optical potential 
of the KVOR model is the closest to that with m^/m7v =0.8 of [17]. Since 
in application to HIC particle spectra should be described in a large energy 
range (up to and above proton momenta ~ 1 GeV/c) in average the KVOR 
description of the data is optimal. 

The isovector part of the optical potential 11^^^ — Uop^ is less constrained by 
the data, cf. [50]. Therefore we do not consider it here. 

The ratios x^b (see Eq. (!3T!) ) are not well fixed experimentally. Different pos- 
sible choices are reviewed in [5T|52] . Within a quark counting model (cf . " case 
I" of Ref. [IS]) one gets Xmb = 1 for non-strange baryons (in our case A's) and 

2 

Xu;A = X^Y. = X^Y.* = XpE = Xpj:* = XpA = 2x^E = 2XpH = " (57) 

in the strange sector, whereas the constituent SU{6) quark model, gives dif- 
ferent values, e.g., Xp= = 1, Xp\ = 0. In this paper we use the quark counting 
values first and then allow for a variation of these parameters. 
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Fig. 1. Energy dependence of the nucleon optical potential for N = Z,V = 0. 
Results for different values of the nucleon effective mass ([8]) m^/niN =0.54, 0.71, 
0.8 and 0.85 calculated in the Walecka model [17] are plotted by dotted, thin con- 
tinuous, dashed and dashed-dotted lines, respectively. The solid line corresponds 
to our calculations with parameter choice ()5ip . Shaded area shows uncertainties in 
extrapolation from finite nuclei to cold nuclear matter |47| . 

The ratios x„h are determined with the help of the relations (compare with [TB]): 

-^bind('^o) = (glN'iT'o/l'nujino) ml])x^h - {rriN - m%{no)) x^h, (58) 

where £'^ind('^o) is the binding energy for the hyperon h. There is a convinc- 
ing evidence from the systematic study of hypernuclei that for A particles 
^bind - -30 MeV. For E we adopt the value Eg^^ ~ -18 MeV and for S 
hyperon -E^ind — MeV (following "case I" of [18]) is taken. There are 
no experimental data for S*. For S* we take the same value as for S. One 
could use other parameter choices within experimental error bars, e.g. follow- 
ing "cases II-IV" of Ref. [IS]. 

We put the coupling g^^, = 0, because lj does not decay in two pions. As 
follows from the 7t~ atomic data, one needs a slight vr" energy shift upward, 
~ 10 ^ 30 MeV for = Z at the saturation density ub = Uq, cf. |53j . 
The value of the pion S-term estimated from scattering data is rather small: 
E^jv — 30-i-45 MeV. Since both these values (energy shift and S^^r) are small, 
we can assume g^^^ = (^o-tt — 0. The p-wave ttNN and vrA^A interactions are 
disregarded within our RMF-based model, as all other p-wave effects. Thus 
for A^ = Z we deal with free pions. 

The pvr-meson coupling gp.,^, which is necessary to describe the pion behavior 
in isotopically asymmetric matter, can be found by matching with the s-wave 
Weinberg- Tomazawa term of 7r~ polarization operator 
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Qp^Rq ^ nn/{Vl) for ^ 0, Z = 0, (59) 

as it is motivated by the chiral symmetry, f^^ = 93 MeV, cf. Thus one 
gets 



9pn ^ . (60) 

The corresponding value gp.„ ~ 6 is consistent with that follows from the 
universality relation, cf. [25] . 

For the kaon coupling constants, we take g^K = Qi^n /"^^ 9pK = \9pN as follows 
from quark counting, cf. [51], and g^K is evaluated from the and K~ 
optical potentials at = no : 

Uopt^^ \no, T = 0) = -g^K cr{no, T = 0) ± g^K ^^oino, T = 0) . 

The experimental values of U^^{no,T = 0) are in the range ~ 20 ^ 30 MeV. 
It is known that the T^kn sigma term is significantly larger than T^t^n. How- 
ever the former quantity is not well determined and may vary in a broad 
range, T,kn ~ 150 -h 400 MeV. The values of a deep potential {no,T = 
0) ~ —100 -. — 200 MeV are derived in Refs. [51] from analysis of kaonic 
atom phenomenology, whereas self-consistent calculations based on a chiral 
Lagrangian [2U] and coupled- channel G matrix theory within meson exchange 
potentials [55] yield U^^; (riQ, T = 0) ~ -50 ^ -80 MeV. In Ref. [18] the effec- 
tive value Tjkn — 150 MeV is extracted using the analysis of the kaon-nucleon 
scattering of [20]. Thereby, we apply a shallow K~ potential {no,T = 
0) ^ -80 MeV. For the K+ potential we take ?7„^t^(no,T = 0) = 21.4 MeV. 
With these potentials we find coupling constants 



(^) : 9*aK = 0-13 gaN, 9^K = ^9i^N- (61) 

The density dependence of the and is even less constrained. Thus 
we also use an another set of kaon couplings obtained with the help of the 
scaling 



(B) : g*K = 0.13 g^N Xa/Xa{no), 9tK = ^9u,n Xuj/xU^o) (62) 
in accordance with the above scaling hypothesis. 

In Fig. [2] the kaon and antikaon dispersion curves are shown versus the baryon 
density for iso-symmetric system at vanishing temperature for both sets of 
couplings. The difference between the kaon (and antikaon) dispersion curves 
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Fig. 2. Kaon and antikaon energies for p = as a function of the baryon density 
for N = Z,T = 0, V = for two sets of couplings: Solid curves are for the (A) 
set without scaling, Eq. ([6T]) . and dashed curves are for the (B) set with scaling, 
Eq. 

for these two parameter choices is tiny up to 2no and becomes significant at 
higher baryon densities. Energies of antikaons ljk- reach zero at — 4.5no 
and ~ 7no for parameter sets {A) and (B), respectively. If a deeper K~ optical 
potential is used, as suggested in [54], one obtains that uk- = at smaller 
density. We note that the points ujk~ = are not the critical points of a 
antikaon condensation. In HIC strangeness is conserved. Then K~ mesons 
can be created only in pairs with mesons at T = 0, if the hyperon Fermi 
seas are not filled. Then the condensation condition for kaons and antikaons 
looks like ujk± (p = 0) ± fistr = 0. In the case of infinitely long-living matter the 
strangeness is not conserved, and we would deal with the kaon condensation 
at ujx- = for T = 0. The consideration of the dense N = Z system at T = 
has only pedagogical interest, however. 

Effect of the ?7-meson on thermodynamic quantities is minor. In our calcula- 
tions of characteristics of HIC we take S^jy — 140 MeV, K^^r = 0.2 fm (see 
Eq. (fT9!) ). These values are twice smaller than the average- weighted values 
S^AT ~ 280 MeV, k^at = 0.4 fm used in Ref. [30]. We make this choice in order 
to simplify the consideration by avoiding a possibility of rj condensation in a 
wide baryon density - temperature range of our interest. Differences of these 
two parametrization in the i] meson mass are shown in Fig. [3J 

In case when there is no information on interactions of a particle with mean 
a-, uj- and p-meson fields, or if it is known that this interaction is rather weak, 
we treat this particle as a free one. As we have mentioned, we consider ?7'(958), 
i^*(892) and (/9(1020) as free particles. We also consider VL and H* as free ones. 
For S* we use the same couplings as for E. 

Though the constructed model describes arbitrary iso-asymmetric systems, in 
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Fig. 3. Effective rj meson mass given by Eq. ()8ip in Appendix A as a function 
of the baryon density for N = Z, T = 0, for two sets of couplings. Solid curve: 
S^AT = 280 MeV, k^at = 0.4 fm. Dash curve: S^at = 140 MeV, k^at = 0.2 fm. 

this paper we further focus on the study of isospin-symmetric nuclear matter, 
N = Z, and disregard small Coulomb effects. In calculations presented below 
we use by default parameters (1321) . (|^ . unless other is specified. 



5 SHMC EoS for T = and N = Z 

At vanishing temperature our model differs from that of KVOR [12] in two 
aspects: (i) we take into account a possibility of occupation of the Fermi seas 
by different baryon species at higher baryon densities; (ii) we incorporate 
a possibility of condensation of the (quasi) Goldstone boson fields, when it 
occurs. 

The left panel of Fig. H] shows the baryon density dependence of the ratio 
of the effective mass to the bare mass for nucleons and u and p excitations, 
m*j^/mN = mPJ"^* /m^ = rnJp'^^* /mp = $ for $ > 0, see Eqs. ([8]) - ([II]), and 
(l44l) . as well as for a excitations, rn^^"^^* /m„, given by Eq. (143|) for T = 0. 

We observe that the effective masses monotonically decrease up to a minimal 
value at the density ub < nmin,B — Skq and then begin to grow. This is 
a consequence of the fact that within our model the masses depend non- 
linearly on the a field and this dependence is determined within self-consistent 
calculations. Due to this feature the SHMC model EoS is getting stiffer with 
increasing baryon density in the range ub < nmin,B and then it becomes softer 
for > nmin,B- Such a high-density behavior could additionally favor the 
deconfinment phase transition at large densities {ub > nmin,B) at T = if it 
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Fig. 4. Baryon density dependence of the effective-to-bare mass ratio at T = 0, 
N = Z. Left panel: for nucleon- uj-p excitations (solid line), and for a excitations 
(dashed line). Right panel: The same as in left panel but for the A(1232) isobar 
(solid line) and the A(1116) hyperon (dashed line). 

had not yet happened at a smaller density. We could chose model parameters 
in such a way that the value of nmin,B would be smaller, in favor of shifting 
the position of the deconfinment transition to smaller density. However it 
would result in the simultaneous decrease of the maximum neutron star mass 
< 1.9 Mq. The latter may come in conflict with experimental data on neutron 
stars, see 

Note that the non-linear density dependence of the nucleon and a effective 
masses resembles the density dependence of the chiral and gluon condensates 
obtained in Refs. [56|57] . 

In the right panel of Fig. H] we show effective masses of the A(1232) isobar 
and the A(1116) hyperon, as representative examples of heavy baryon species. 
Their baryon density dependence is similar to that for the nucleon, with the 
same value of nmin,B- However the value of ml{nmin,B) /''^h ip 7^ N) is higher 
than m*j^{nmin,B)/mN. 

In Fig. [5] we show the density- dependent total energy per baryon for SHMC 
EoS in comparison with that for the Urbana-Argone (Al8-f5u+UIX*) EoS 
[39] (in the HHJ version of Ref. [H]) and for the ideal gas (IG) EoS. Within 
the IG model we include the same particle species, as in the SHMC model, 
but in this case all mean fields and thus all particle interactions are switched 
off. Thus, in the IG model at T = only nucleon Fermi seas contribute. 
It is seen that difference between SHMC EoS and IG EoS grows strongly 
with the density increase indicating to an important contribution of particle 
interactions. The threshold density for the appearance of the A(1232) isobars 
is ~ 12 no and their Fermi sea again melts for > 20 no (these densities 



25 



1400 




Fig. 5. Baryon density dependence of the total energy per baryon (solid line) for 
SHMC EoS, for Al8+(5u+UIX* in the HHJ parametrization [JT] (dotted line), and 
for the IG EoS (dashed line); T = 0, N = Z . Arrows show the threshold densities 
for the appearance and melting of the A isobar Fermi sea. Dashed-dotted curve 
is a partial contribution of the repulsive vector field term of the SHMC EoS. For 
comparison the kinetic energy En /nB — for the SHMC model is shown by the 
thin continuous line. 




Fig. 6. Baryon density dependence of pressure for symmetric nuclear matter at T = 
for the SHMC (solid line) and IG (dashed) EoS. The shaded area corresponds to 
the experimental constraint derived from HIC |59j . 
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are shown by arrows in Fig. [5]). Appearance of A's almost does not affect the 
total energy per baryon. Hyperons do not occur at all for = Z in contrast 
with /^-equilibrium matter, cf. [TB] . 

Though baryon and meson masses begin to increase for Ub > nmin,Bi this 
only moderately affects the stiffness of the EoS, since the suppression of the 
nucleon kinetic term (thin solid line in Fig. [5]) is largely compensated by the 
increase of the repulsive vector meson term (dashed-dotted line in Fig. [5]). Due 
to that SHMC EoS remains stiffer for n > rimin^B compared to the HHJ EoS 
(dotted line in Fig. \5^. 

The SHMC EoS begins to differ from the HHJ EoS for > 4no and this 
difference increases with increase of the baryon density. Such a behavior, cf. 
|12] . results in an increase of the value of the maximum neutron star mass, 
Mmax = 2M0, that is in agreement with the value M^ax = (2.1 ± 0.2)Mq 
(at the la confidence level) derived in [58j from the observations of the PSR 
J0751-I-1807, a millisecond pulsar in a binary system with a helium white 
dwarf secondary. 

In Fig. E] the pressure calculated in the SHMC model (solid line) is compared 
with the experimental constraints coming from the analysis of elliptic flow in 
HlC [59]. As it has been argued in only the EoS with pressure curves, being 
close to the upper boundary of the band, satisfies the maximum neutron star 
mass constraint. Pressure within the IG model of EoS (dashed line) does not 
fulfill the HIC flow constraint. To satisfy the flow constraint at T = 0, one 
definitely needs a much stiffer EoS than that given by the IG model. 

Transport calculations [60] have demonstrated that subthreshold produc- 
tion may provide an important information to constrain the EoS of the warm 
symmetric nuclear matter for ~ 3 uq. Within the last decade the KAOS 
collaboration at GSI performed measurements of the kaon production |61j . 
Analysis of the data [52] led to a conclusion that the EoS satisfying the kaon 
data is compatible with the above-required flow constraint. Both constraints 
hold true with the EoS of the Urbana-Argonne group (A18 + 6v + UIX*) and 
with the KVOR-based SHMC EoS used here. 



6 SHMC EoS for T ^ and N = Z 

6.1 Density-temperature dependence of effective masses of excitations 

The effective masses of the nucleon and uo/ p excitations follow the same scal- 
ing law and coincide, see Eqs. (ITT]) . (Hljl . As it is seen from the left panel 
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of FiglTl at the baryon density ub < Q ^ 8 uq and T < 190 MeV, the effec- 
tive masses of the nucleon and u/p excitations and the a excitation decrease 
when the density grows up and then they start to increase at higher densities 
similarly to the case T = 0. 

Generally, different phase states can be realized within the SHMC model. At 
some density ub = n]^ and temperature T = T'^'^ the a excitation mass may 
reach the value 2m„. Then the decay cr — > 27r becomes forbidden at higher T 
and fiB- Ref. [63] argued that due to long-scale field fluctuations the scattering 
length of two pions at rest should go to infinity at ra^ ^ n'^ , identically to 
the so called "Feshbach resonance" at zero energy to be used in atomic physics 
for cold trapped atoms. The result is also known as a new "strong coupling" 
regime of matter which manifests a liquid-like behavior [M]. However these 
interesting questions are beyond the scope of the present paper since in the 
SHMC model the particle excitations have no widths. Therefore we continue 
to apply our model without any modifications also for Ub > n"^ . 

In addition, the a excitation mass reaches the value at some density n^''' 
and temperature T = T'^^"'. Note that in models with inherent chiral symmetry 
the a and vr masses meet at the chiral symmetry restoration point. As it is 
seen from Fig. [TJ with the SHMC model for the T = case the a excitation 
mass remains always higher than the double pion mass. But at sufficiently high 
temperature, T > 190 MeV, both points, ub = n"^ and n^'"', are reached. 




^^^0 1, MeV 

Fig. 7. The effective masses of the nucleon, p excitations (thick hnes), and of 
a meson excitations (thin lines), N = Z. Two dotted straight hnes show the a 
mass levels rrio- = 2771,^ and rrio- = ttItt (bottom). Left paneh The baryon density 
dependence at different temperatures. Right paneh Temperature dependence at the 
baryon density = (continuous hnes) and ub = 5no (dashed hnes). 

As follows from the left panel of Fig. [71 for the phase range T < 190 MeV and 
nB < (6 -i- 8) uq the density behavior of effective masses is similar to that in 
models with inherent partial restoration of the chiral symmetry. The values 
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of effective masses decrease at low densities with the temperature increase in 
accordance with behavior of $ = 1 — /. For T > 190 MeV all the effective 
masses are getting small rather sharply but then grow slowly with the increas- 
ing density. At temperature T =188 MeV the critical value n'^ is reached first 
at n'^ ~ Auq and this phase state is left at n'^ ^ Sns (see crossing of the 
thin dashed line with the horizontal one). At T = 200 MeV the mass of the a 
field is less than ttItt for < < IOuq. 

As it is seen from the right panel of Fig. [3, the temperature dependence of the 
effective masses of nucleon and sigma excitations is small up to T ~ 170 ^ 180 
MeV. For higher temperatures the effective masses begin to decrease abruptly 
and for ub = we have T*^^ ^ j^cha ^ j^gg ]\/[gY jf proceeds to the 
dense matter (n^ = 5no) the difference between these temperatures is about 
few MeV. Within this narrow temperature interval the second derivative of 
the effective mass with respect to temperature changes the sign. The effective 
nucleon-cj-p and a excitation masses reach zero at the same critical temper- 
ature Tea, about 210 MeV. Since the coupling scaling functions Xo- and Xuj 
follow the same dropping trend as the mass scaling function $, in vicinity of 
Tea we deal with a gas of almost massless excitations. Similar result has been 
obtained in [53] using a generalized local symmetry approach and vector man- 
ifestation arguments. The anti-nucleon yield rapidly increases due to a sharp 
decrease of the nucleon mass. 

For T > Tea the effective nucleon mass becomes negative. For the first time 
such a behavior of the nucleon effective mass has been found within the stan- 
dard RMF model including A resonance in [66]. Authors suggested a specific 
choice for the resonance-cx couplings that allows to restore positiveness of 
masses. Actually in the region where the effective nucleon mass is negative 
nothing dramatic happens. The nucleon spectrum given by Eqs. (127|) . fl28l) 
continues to be well defined since these equations enters rather than m|^r. 

Note that in our model the effective u and p excitation masses follow the law 
(fTTl) . (l44l) . Thus they only touch zero at T = Tea and become again positive 
for T > Tea- Therefore u and p condensates should not appear at T > Tea- 

However the effective mass of the a excitation {m^'^^^*, see Eq. fl^ ) in Ap- 
pendix A) is getting imaginary at T > Tea- Thus, the ground state proves to 
be unstable with respect to the Bose condensation of the a excitation field. 
The stability is achieved due to the self-interaction between the cr-particle 
excitations, see Appendix B. Such a condensation might be called "hot Bose 
condensation" since it occurs at T > Tea, in contrast with the standard Bose- 
Einstein condensation appearing with the temperature decrease. A possibility 
of hot Bose condensation has been considered in [23], within a different model 
which includes effects of particle widths. In order not to complicate consid- 
eration we avoid description of temperature region above Tea in the present 
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work. 
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Fig. 8. Baryon density dependence of energies ujk (upper panels)) and ijJk — i^s 
(lower panels) for kaons (thick lines) and antikaons (thin lines) at p = at three 
values of the temperature T =42, 132 and 188 MeV, N = Z. Left panel: the {A) 
set for kaon couplings without scaling, Eq. ()6ip . Right panel: for the {B) set of 
couplings with scaling, Eq. ([62|) . 



Fig. [8] (upper panels) presents the baryon density dependence of the energy 
for kaons [K^ or K'^) and antikaons {K~ or K^) with zero momenta at three 
values of temperature, see Eqs. (l77 l) -( 179l) of Appendix A. Kaons {K^ and 
K^) (as well as antikaons {K~ and K^)) have the same dispersion relations 
in the iso-symmetrical matter, if a small Coulomb contribution is neglected. 
We see that the energy only slightly depends on the temperature for 
T < 170 MeV, ujk+ decreases with the T increase. Temperature dependence 
of the K~ energy is also minor for T < 170 MeV, ujk- increases with the T 
increase. Thus for T < 170 MeV the density dependence remains similar to 
that at T = 0, see Fig. [21 For higher temperatures the T dependence becomes 
significant for both kaons and antikaons. At = 0, branches ujx+ and ujk- 
coincide. For set B of scaled couplings the dispersion curves are more flat than 
for set A. For set A, ux- vanishes at ~ (3 -i- 5)?t,o (for T < 188 MeV) and 
for set B, for ~ 7no. However it does not mean that the condensation 
occurs. The necessary condition for condensation is ujk{p = 0) — /i^ = 0, 
fJ's = Z^str for kaons and fig = — /istr for antikaons, see Eq. (!83|) of Appendix 
A. This difference is plotted in the bottom part of Fig. [HI It is seen that the 
kaon condensation condition is never fulfilled. Antikaon condensation takes 
place only if the density-dependent scaling is neglected (set A), at ub ~ Suq 
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and T ~ 190 MeV. With set B of couplings the antikaon condensation does 
not occur in the relevant density-temperature range. Therefore we perform 
subsequent calculations of HIC for the case B. 

1,0 
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E 0,4 
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2 4 6 8 10 12 14 

Fig. 9. Baryon density dependence of the r] effective mass at four values of the 
temperature T =42, 132, 188 and 200 MeV, N = Z, for S^at = 140 MeV and 
l^-rjN = 0.2 fm. 

Fig.[9]shows that the rj effective mass monotonically falls down with the baryon 
density if the temperature is not very high. However at T ~ 200 MeV the rj 
effective mass becomes almost independent of ub for < 12no, m*/mr, — 0.4, 
within the ub — T region depicted in Fig. [91 For > 15no the in-medium 
rj mass practically is independent of temperature. The rj condensation does 
not appear for T < Tea- If we used parameter choice S^tv = 280 MeV and 
fir)N = 0.4 fm (see solid line in Fig. [3]) we would meet with the t] condensation 
problem. 

6.2 SHMC EoS for baryonless matter. 

Now let us consider the case — that is close to conditions realized at 
RHIC. The decrease of the hadron masses with increase of the temperature 
for ^0 has been found in [21] as the consequence of the blurring of the 
baryon and meson vacuum. Here we obtain a similar effect but within the 
quasiparticle picture. 

In Fig. [TU] temperature dependence of the pressure is shown in units (left) 
and the specific heat, in units (right) for = 0. On the left panel the 
solid curve presents our calculation with parameters determined in sect. HI 
The dashed curve shows pressure in the case when contribution of all baryons, 
except neutrons and protons, is artificially suppressed while the dotted curve 
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T, MeV T, MeV 

Fig. 10. Temperature dependence of the reduced pressure (left) and reduced specific 
heat (right) for iso-symmetric baryonless matter. The solid line presents our full 
calculation. The dashed line shows the case when contribution of all baryons besides 
neutrons and protons is artificially suppressed. Two dashed-dotted lines are the case 
when all gmb couplings except for nucleons are suppressed by factors of 2/3 and 1/3, 
respectively. Squares show the lattice QCD result for the 2+1 flavor case [67j. The 
dotted curve is for the IG model. 

is given for the IG model. All curves are close to each other for T < 160 + 
170 MeV. Differences start at higher temperatures where high-lying baryon 
resonances come into play. At T = Tea the effective mass squared for the a 
excitation drops to zero and becomes negative for T > Tea- In order not to 
complicate the consideration we avoid discussion of the hot Bose condensation 
of cr excitations. Therefore we cut the solid line at T > Tea- In the case when 
baryon resonance contributions (except nucleons) are artificially suppressed, 
the effective nucleon mass as well as the a excitation mass approach zero only 
at T — i> oo (if ?7 is a free particle, otherwise rj may condense at Tcr^ — 265 MeV). 

We also compare our result with the lattice QCD calculations for 2 + 1 fla- 
vor [67]. One usually believes the lattice data for the quark-gluon sector 
(i.e. for T > T^ec ~ 175 MeV and for ub — 0). On the other hand, there 
are doubts that the lattice calculations produce appropriate results for the 
hadronic sector, T < T^cc, since they get unrealistically high value for the 
pion mass, ~ 700 MeV, instead of the physical value 140 MeV, cf. |67j . 
Recently a higher deconfinement temperature was obtained in lattice QCD 
calculations for similar system with almost physical value of the pion mass, 
^dcc = 192(7)(4) [68]. However thermodynamical characteristics still have not 
been recalculated. Therefore we may use any hadron EoS for T < Tdcc (with a 
not precisely known value of T^cc) without referring to the lattice data. For a 
larger temperature, one could expect that the quark phase becomes energeti- 
cally preferable. It happens if the pressure in the hadron phase is less than in 
the quark phase. Oppositely, Fig. [in] shows that the pressure extracted from 
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the lattice calculations is significantly smaller than that obtained within the 
SHMC model for the hadron phase at T > 170 MeV. This may mean that 
we have no deconfinement phase transition with our EoS with the values of 
the coupling constants used. Instead, within our model we obtain a state of a 
high-temperature hadron gas of many baryon resonances (quasiparticles in the 
given model) and bosons, with small effective masses. This phase is enriched 
by antiparticles. As we have mentioned above, this result coincides with that 
of Ref. [23], though here it is obtained within a phenomenological quasiparticle 
scheme while in [23] it is a consequence of blurring of the hadron vacuum. 

The deconfinement transition could be constructed if the (lattice) quark-gluon 
EoS is matched with our EoS for lower T (Tdcc ^160 MeV) when our hadron 
phase pressure is not yet too high (compare solid and dashed lines in Fig. [TOj) . 
Another possibility for the deconfinment transition can be associated with 
a different mechanism: The overlapping of the hadron cores, if hadrons are 
considered as composite particles. We have a dramatic increase of hadron 
degrees of freedom at T > 170 MeV. Thus the hadron cores may become 
overlapping for such temperatures. If this mechanism works, the deconfinment 
transition should be treated as an enforced Mott-like transition occurring due 
to the melting of composite hadrons, rather than matching the pressures of 
two different phases. 

In Fig. [To] solid lines present results of calculations with the default param- 
eters defined in sect. 4. Here we have Tea — 210 MeV. Dashed-dotted lines 
demonstrate the cases when gmb couplings, except for nucleons, are suppressed 
by factors 2/3 and 1/3 (as indicated on the plot). The latter case (with 1/3 
prefactor) allows to fit the lattice data up to T ~ 240 MeV. In this case a 
quark liquid would masquerade as a hadron one. In principle one could fit 
the lattice data in a still larger region of temperatures (e.g., up to 500 MeV) 
introducing x^- < $. A violation of the universality of the a scaling would be 
in a line with that we have used for u and p, 77^^ 7^ 1 and ?7p 7^ 1. However 
we will not elaborate this possibility in the present work. The dashed-dotted 
curve labelled by 1/3 is cut at Tc^ ~ 265 MeV. At this point the ?7-meson mass 
becomes imaginary within our parametrization and rj condensate arises. If 77 
were treated as a free hadron, we would obtain T^a — 330 MeV in this case. 

The right panel of Fig. [10] demonstrates the temperature behavior of the spe- 
cific heat for ub = 0- With the standard choice of couplings (solid curve) 
we observe a sharp peak with maximum at T ~ 190 MeV. At this point the 
second derivative of the effective nucleon-tu-p and cr-excitation masses changes 
the sign. The specific heat retains a continues function. The obtained behav- 
ior is typical for the strong crossover transition. Note that such a behavior of 
the specific heat was also found within the standard RMF model [69]. In the 
case of the second-order phase transition, the specific heat would be discon- 
tinues at the critical point. When gmb couplings are suppressed, the peak is 
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smoothed that reminds about a weak crossover. Often the temperature at the 
maximum of the specific heat is associated with the critical temperature of a 
phase transition. However note that in our case the position of this maximum 
T ^ 190 MeV shghtly differs from the T = T^^^^ and it differs also from Tea 
(~ 210 MeV), see Fig. [71 This fact of non-unique value of the critical tem- 
perature Tc is also manifested in the lattice calculations: Analysis of different 
thermodynamic quantities leads to different numerical values of Tc even in 
the continuum and thermodynamic limits [70]. One should keep in mind that 
there is no liberation of internal (quark-gluon) degrees of freedom of hadrons 
in the RMF models. Note that the P/T^ curve calculated in the IG model is 
below our result even if only nucleons are taken into account (compare dotted 
and dashed curves in left panel of Fig. [TOl) . This is due to the decrease of the 
nucleon mass with increasing temperature in the SHMC model. The specific 
heat Cv/T^ in the IG model does not saturate for high T, whereas all curves 
of the SHMC model tend to constant values. 



6.3 SHMC EoS for hb ^ 



In Fig. [TT] we show temperature dependence of the pressure (left) and spe- 
cific heat (right) for baryon densities n^/no =2, 5 and 10. The pressure gets 
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Fig. 11. Temperature dependence of the pressure (left) and specific heat (right) at 
ub ^ for N = Z for three values of . 

minimum at temperature T ~ (175 ^ 200) MeV and grows for higher T. The 
calculations are stopped at T = T^a as was mentioned above. The value Tea 
depends on only moderately. We see that the peak of the specific heat 
survives at finite baryon density. The value of the temperature T^ax corre- 
sponding to the peak is slightly shifted up with the baryon density but the 
hight of the peak changes significantly. In particular, the position of the c^/T^ 
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maximum is Tmax ~ 190 MeV at ub ^ 2no and moves to about 200 MeV for 
the density ub = lOriQ. 

6.4 Particle densities 

The particle density of species i is given by 



where the spin- degeneracy factor Qi = (2sj + l) and is the particle occupation 
function, see Eqs. (l28D-([iD,([82]),([MD. In the SHMC model the spectra of most 
particle species are getting softer with the baryon density increase because of 
in- medium effect. This occurs for all baryons and for u and p due to the scaling 
of the effective masses, as well as for a and (quasi) Goldstone excitations {K~ , 
and 77) as a consequence of their interaction with mean fields. Therefore 
the densities of these particle species are larger than in the IG case. 

Temperature dependence of the density for various particle species (and for 
their antiparticles) for iso-symmetric matter at different values of baryon den- 
sities is shown in Fig. [T^i In the = case the particle and antiparticle 
yields coincide. As it is seen, all the species besides nucleons exhibit very simi- 
lar behavior (compare solid and dashed lines). Particle number density slowly 
grows with the temperature increase till T ~ 170 MeV and then rapidly goes 
up. In the range T < 170 MeV the SHMC results are rather close to those of 
the IG model (except for a high density, see example hb = Srag) but drastically 
diverge at higher temperatures. This difference is naturally explained by the 
rapid decrease of in-medium masses in the SHMC model at these temperatures 
(see Figs. H|7|) . In contrast, at T < 170 MeV the nucleon density un stays 
almost independent of temperature aX ub — O.Srio, or significantly decreases 
at higher ub- These facts are due to a stabilization effect of the baryon con- 
servation law and a strong growth of production of baryon excited states and 
hyperons with the baryon density increase (cf. A and A particle densities in 
Fig. [12]). Similarly to other species, the nucleon particle density rapidly goes 
up with further temperature increase (for T > 170 MeV). 

Temperature-density behavior for antiparticles is quite different. In spite of 
the fact that the iso-symmetric system is considered, even for IG the densities 
of created kaons and antikaons are different (except for the = case). 
This is a consequence of the total strangeness conservation which takes into 
account also hyperon species. Thus /igtr proves to be non-zero even for N = Z 
that results in nK+ 7^ nx-- Although the correction to the energy dispersion 
curve is positive (repulsion) for kaons and negative (attraction, see Fig. [2]) 




(63) 
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Fig. 12. Particle number density in no units versus temperature at different baryon 
densities for A, N, A and their antiparticles, N = Z. Our model results are 
plotted by solid line for particles and by dashed-dotted line for antiparticles. Dashed 
and dashed-double-dotted curves correspond to particles and antiparticles in the IG 
model. 



for antikaons, their number densities intersect at temperature 160 J; T < 
180 MeV for 7^ and then, with the subsequent temperature increase, 
the antikaon density even exceeds that for kaons. But it is only apparent 
violation of strangeness conservation since excited kaon states K* contribute. 
Because K* is treated as a free particle, nK*+ ^ uk*- at high temperature 
and baryon density. In particularly, at T = 190 MeV and ns = Srio we have 
hk = 0.46no, riji = 0.68?7,o and for the excited kaon state uk* = l.lSno, 
n/f* = 0.07no. Indeed, nx < but the total strangeness of these four kaon 
species is ~ O.Suq which is compensated by hyperons. As for antinucleons, 
antideltas and antihyperons, their behavior looks very similar: the yield of all 
antibaryons is markedly suppressed at T < 170 MeV and abruptly grows up 
at higher temperature, T > 180 MeV. This increase is not reproduced by the 
IG model which produces significantly lower hadron densities. Generally, the 
baryon density dependence of the particular number density seems to be not 
as strong as temperature one. 
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7 Application of the model to HIC 



The SHMC model describes the EoS of hot and dense hadronic matter in a 
broad range of temperatures and baryon densities. In HIC the dense mat- 
ter created in the initial stage is expected to rapidly thermalize and then 
it expands without significant generation of the entropy S. In the process of 
expansion some particles may leave the fireball carrying away a part of the en- 
tropy. Thus more appropriate characteristic which is approximately conserved 
is the entropy per participant nucleon (S/Nb). This thermodynamic quantity 
should be conserved in quasi-equilibrium case and it is also less affected by 
any possible particle loss or gain from the fireball during the expansion stage. 
The predictions of our models for the evolution path in the (T, yU,B)-plane as 
obtained from the SHMC EoS under condition of the fixed S/Nb are shown 
in Fig. M 
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Fig. 13. Isentropic trajectories for central Au+Au collisions at different bombarding 
energies calculated in our model (solid lines) and IG (dashed lines). Experimental 
points with error bars are taken from [73]. The freeze-out points marked by stars are 
obtained in this work (see the text below). Thin line corresponds to the freeze-out 
curve in [73] while dash-dotted line is that from |72j . 

For the initial fireball expansion stage, in the case of semi-central Au+Au 
collisions at different bombarding energies below the top AGS energy, the re- 
duced entropy ratios, S/Nb, were estimated within a transport Quark-Gluon 
String Model (QGSM), as described in |34p1 . The temporal S/Nb dependence 

^ As shown recently [71] dynamical trajectories described by the QGSM are quite 
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exhibits some saturation [71] , values being taken as an input for our isentropic 
calculations and indicated in Fig. [131 The presented isentropic curves follow 
exactly and unambiguously from the EoS. Some uncertainties in location of 
these curves are coming from uncertainties in the estimate of the reduced en- 
tropy related to the bombarding energy. In our case 6{S/Nb) ~ 0.1 and it 
unessentially influences the trajectory location. The {T,fiB) trajectories cal- 
culated within the SHMC model show turning points those positions correlate 
roughly with the freeze-out curve. This fact was noticed earlier in |7l] where 
the EoS with a phase transition was used. There the high-temperature part 
of trajectory with dT/dfiB > was associated with the quark-gluon sector of 
EoS but in the present work it is due to a strong decrease of hadron masses in 
the considered area of the phase diagram. Note that there is no turning point 
in the IG case. 

Two chemical freeze-out curves are shown in Fig. [T31 Dashed-dotted curve 
corresponds to the IG EoS with the chemical freeze-out condition that the 
energy per hadron equals to 1 GeV {< E > / < N >= 1 GeV) [72]. It is of 
interest that in the considered {T,^b) range this freeze-out curve is also very 
close to that obtained for the IG EoS with the fixed baryon density ub = 0.12 
fm~^ [73]. So, the net baryon density of states above the dashed-dotted line 
in Fig. [13] is ub > 0.8no. Maximal densities explored in this phase diagram 
range from ~ 2.3no to about 4no when the bombarding energy increases from 
2 AGeV till the top AGS energy (these numbers depend on the region over 
which averaging was made). The thin line is obtained by interpolation of 
the (T, fis) fitting parameters extracted at every available bombarding energy 
by the minimization of the difference between experimental and theoretical 
(calculated within a statistical model with the IG EoS) hadron abundance |73j . 
In any analysis the thermodynamic quantities (T, /is) at the freeze-out are de- 
rived from the analysis of measured particle ratios. The straightforward con- 
sequence of the statistical assumption is that the mean hadron multiplicities 
should be calculated in the full (47r) phase space. However, due to limited ex- 
perimental acceptance it is not possible to do that in each case and instead the 
particle ratios at the middle rapidity are used. Making use of middle rapidity 
ratios implies that the measured distribution dN/ dy is approximately constant 
over the same range. The analysis of the expected dispersion of rapidity distri- 
butions shows that the use of full space multiplicities is better suited over the 
energy ranges of AGS and SPS [75]. Unfortunately in the energy range con- 
sidered in Fig. [13] the Att particle ratio measurements are available only at two 
energies, and the thin line is obtained by the (T, /x^) interpolation based on 
the analysis of the middle rapidity data. Generally, both heuristic freeze-out 
curves describe quite well the extracted (T, /i^) values and they differ more 
noticeably in the presented range Eiab ~ 10 AGeV, below the top of AGS en- 



close to those in the UrQMD model providing the closeness of initial fireball stages 
in these two approaches. 
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ergy. Before to explain how freeze-out points marked by stars in Fig. [13] were 
obtained in the given work, we consider which hadron abundance is predicted 
by the SHMC model and which values of thermodynamic parameters may be 
inferred from comparison with experiment. 
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Fig. 14. Particle ratios calculated along the isentropic trajectories for central 
Au+Au collisions at Eiab = 4 AGeV. Left panel: Solid and dashed lines are cal- 
culated right before freeze-out for the SHMC and IG EoS, respectively. Right panel: 
Solid line is calculated in the SHMC model right after freeze-out has occurred, 
with taking into account the canonical strangeness suppression. Dashed line is that 
without strangeness suppression. Shaded bands correspond to uncertainties of ex- 
perimental data [73]. The dotted vertical line is our estimate of the freeze-out tem- 
perature. 

In Fig.[T4]the Eiat = 4 AGeV case is exemplified in detail. In the left panel the 
particle ratios for the most abundant hadrons are calculated along isentropic 
trajectories of the SHMC and IG EoS's depicted in Fig. [131 Note that these 
quantities are in-medium ratios which would be "inside" hot and dense matter 
right before the freeze-out, but measured ratios correspond to free particles 
after the freeze-out. Although trajectories are remarkably different for the IG 
and SHMC EoS's, the particle ratios turn out to be not so sensitive to the 
choice of EoS (besides the K~ /ti~ ratio), since freeze-out baryon densities are 
rather low. In the right panel of Fig. [H] the particle ratios are shown just after 
the freeze-out. 

Up to now there is no appropriate theory of the freeze-out though there exist 
many different recipes. Generally, transition from the collective expansion to 
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kinetic stage and then to free particle streaming should be continuous and it 
takes some finite time. For the sake of simplicity some sudden approximation 
is often applied. If we assume a prompt freeze-out concept, the observable 
yield will be defined by in-medium particle spectra dN/(fp, as we calculated 
them within the SHMC model, but additionally multiplied by a pre-factor 
^/^^p"qr^/ y/m*'^ _l_ p2 quasiparticle undressing, cf . [76] . However the 

assumption of prompt freeze-out might not be applicable at least for some 
particle species. Besides, one should take a especial care about the total energy 
conservation. 

In Ref. [77] another more conventional choice was suggested. At crossing the 
freeze-out point (hypersurface in a general case) the change from the in- 
medium to IG EoS occurs in a "shock-like" way: One demands the energy, 
momentum, and the net baryon charge and strangeness conservations. We 
do not consider the fireball expansion dynamics but use the isentropic tra- 
jectory with our EoS. Thus we know nothing about momentum conservation 
and in our simplified case only conservation of energy and charges is taken 
into account. Certainly, important collective fiow effect is out of considera- 
tion. Therefore below we focus on analysis of particle ratios, where this effect 
is cancelled, cf . [721175] . We use this scenario [77] in our calculations presented 
in Figs. [14] (right panel) and Fig. [151 In this case the attractive in-medium 
interaction brings to an increase of temperature of free gas after freeze-out. 
As shown in the right panel of Fig. [M] (dashed lines), such a procedure results 
in the agreement with experiment of particle ratios at temperature by about 
10 MeV higher than that inside the medium (cp. the left panel of Fig. [T^ . 
Certainly, for the IG EoS there is no changes due to the energy-momentum 
conservation. 

All above-mentioned calculations have been done in the grand canonical en- 
semble. However for description of strangeness production at not too high 
temperatures when a number of strange particles is small, this approach is 
not quite appropriate and the canonical ensemble for strangeness should be 
used [78]. Replacement of the grand canonical description by the strangeness 
canonical one results in extra temperature-dependent suppression factor for 
strange particle densities which reaches unity for T > 100 MeV. The factor 
was calculated in the standard way [75][75] . This canonical strangeness sup- 
pression effect is clearly seen from comparison of solid and dashed lines in 
Fig. [m (right panel). Note that, if both effects are taken into account, the 
best agreement of calculated particle ratios simultaneously with all measured 
ones is reached at the temperature T ~ 78 MeV which can be treated as a 
freeze-out temperature Tf^ in the given case (shown by the vertical line). 

Final SHMC results for hadron ratios at other four AGS energies are presented 
in Fig. [151 Calculated along isentropic trajectories, these ratios take into ac- 
count both the shock-like freeze-out and canonical suppression effect (compare 



40 



140 




T, MeV 



Fig. 15. Particle ratios for central Au+Au collisions calculated at different collision 
energies. Solid lines are SMHC results at the shock-like freeze-out with taking into 
account canonical strangeness suppression. Shaded bands correspond to uncertain- 
ties of experimental data [73]. In the case of Eiab = 2 AGeV the experimental line 
for K~ /'K~ shows only the upper limit of this ratio. The dotted vertical lines are 
freeze-out temperatures. 



with in-medium ratios shown in Fig. [H] for the case Eiab = 4 AGeV). One can 
see that at every bombarding energy it is possible to fix some freeze-out tem- 
perature Tfr by condition of the best agreement of the calculated ratios with 
experiment. The extracted Tfr are shown in Figs. [HI [15] by vertical dashed 
lines. One should note that the proton yield includes both direct protons and 
those feeding from the resonance decays. However in the intermediate energy 
range considered, a part of protons is bound into light complex particles (rf, a). 
This nucleon coalescence effect is the higher, the lower the bombarding energy 
is, and this effect is disregarded in our model. Thus a discrepancy with exper- 
iment for p/yr^-ratios seen in upper panels of Fig. [T^ should not be taken too 
seriously. 

If the Tfr on the isentropic curve is known, an appropriate baryon chemical 
potential for the IG EoS at the freeze-out, /x/^, can be found. These pairs of 
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{Tfr,fifr) are indicated by stars in Fig. [131 The new freeze-out points, based 
on the same middle rapidity particle ratios, slightly differ by some shift in 
the chemical potential from those obtained by direct fitting of these ratios 
in statistical theory [73]. As to Tfr, our freeze-out temperature for Eiat = 10 
AGeV is noticeably below the fitted one. In contrast with lower energies, at 
the top AGS energy the set of 10 particle ratios was measured and used in 
the statistical model fit, while our results are based on 4 ratios, as shown 
in Fig. [151 If the d/p,p/p,A/A,ip/K~^ ratios are excluded from this analysis, 
then one gets Tfr = 108 ± 9 MeV and hb = 555 ± 18 MeV [73] what is in 
a reasonable agreement with our result (see stars in Fig. [T31) . Note that this 
new method for deriving freeze-out points in the phase diagram allows one to 
keep some memory on the collision dynamics (the value of the entropy per 
baryon and isentropic trajectories inherent to the final expansion stage) and 
takes into consideration the in-medium particle modification. The freeze-out 
point is taken on the phase trajectory which depends on the EoS as is seen 
from comparison between the SHMC and IG models in Fig. [131 Note that 
partial ratios are sensitive to the choice of the freeze-out scenario. If we used 
the prompt freeze-out concept [76] we would obtain other yields. 
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Fig. 16. Isentropic trajectories for ultrarelativistic central Au+Au collisions at 
different bombarding energies calculated in the SHMC model with the default set 
of parameters (solid lines) and with the suppressed couplings by the factor 1/3 
for all baryons besides nucleons (dashed lines). Filled diamonds and triangles are 
obtained from the 47r particle ratios in [73] and [75], respectively. Filled circles are 
RHIC data based on the middle rapidity particle ratios [73]. Open circles, squares 
and triangles are the lattice 2-flavor QCD results [81] for S/Nb = 30, 45 and 300, 
respectively. Two freeze-out curves are the same as in Fig. [T3l 
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Let us extend our analysis to higher bombarding energies. The entropy per 
baryon participants was calculated in [80] within the 3-fluid hydrodynamic 
model assuming occurrence of the first order phase transition to a quark-gluon 
plasma. The energy range from AGS to SPS was covered there. We use the 
S/Nb values for Eiab =158, 80, 40 and 20 AGeV, at which the particle ratios 
were measured by the NA49 collaboration (cf. [73j). At the RHIC energy 
we put S/Nb = 300 in accordance with the estimate in [81]. Note that in 
this analysis the {T,^b) thermal parameters were determined using the 47r 
particle ratios, besides the RHIC point. The difference between two sets [751175] 
of (T, /ifi) freeze-out points is caused by more elaborated statistical model 
implemented in Ref . [75] . It is of interest that both sets strongly correlate with 
the freeze-out curve < E > / < N >= 1 GeV [72] and markedly differ from 
interpolation of the freeze-out points based on the middle rapidity particle 
ratios. 

The calculated trajectories for isentropic expansion are presented in Fig. [TBI 
The trajectories (solid lines) behave very similar to those at lower bombard- 
ing energies (shown in Fig. [T3|l exhibiting some flattening above the freeze-out 
curve near the expected phase boundary. Such a behavior reminds the one 
obtained in [S2] where a smeared first-order phase transition is assumed. If 
baryon couplings (besides nucleons) are suppressed by a factor 1/3, the cal- 
culated trajectories turn out to be very close to the lattice QCD results [81j . 
Note that dealing with the finite chemical potential we use here the same 
choice of parameters gmb as in the case ub = (see Fig. [TOj) . 

In principle, we could repeat our particle-ratio analysis in the SPS energy do- 
main for both parameterizations of the SHMC model, with suppressed and not 
suppressed couplings gmb- However for this aim the used SHMC model basis of 
hadron species should noticeably be enlarged by inclusion of higher resonances 
since available set of experimental data for particle ratios is significantly larger 
at high energies [73j . 



8 Conclusions and perspectives 

In this paper the modified relativistic mean-field a-u-p model with scaled 
hadron masses and couplings (SHMC model) formulated for T = in Ref. [12] 
is generalized to finite temperatures. Besides nucleon and mean fields the 
model includes low-lying baryon resonances and their antiparticles, boson ex- 
citations (following the SU{3) concept) and a- u- p-excitations on the ground 
of mean fields. The EoS for T = satisfies general constraints known from 
atomic nuclei, neutron stars and those coming from the flow analysis of HIC 
data. Like in the KVOR model [1], we assume that a, u, p field mass terms, 
as well as nucleon masses, decrease with increase of a combination of a a 
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mean field / = QaNXi.^)^ I'^n (see Eq. (fTOl) ) corresponding to the change of 
the chiral condensate density. The model supposes the simplest choice of the 
Brown-Rho scaling when the change of all masses mentioned above follows the 
same universal law. In order to describe properly the EoS, a similar scaling of 
the coupling constants is introduced (with a slight violation of the universality 
of the scaling law). 

It was shown that at T = our model simulates the (partial) chiral sym- 
metry restoration with the baryon density increase. The baryon-cu-p and a 
excitation masses fall down with the baryon density increase in the interval 
Q < ub < nmin,B — Suq. To Certain extent, this is naturally predetermined by 
the imposed Brown-Rho scaling, in terms of the factor $ = 1 — / (see Eqs. ([8]) 
- dUD)- For higher densities the masses begin to grow up. Such a behavior 
resembles the density- dependence of chiral and gluon condensates obtained in 
Refs. ^m7\ . 

Although the Lagrangian of the SHMC model does not respect chiral sym- 
metry the model simulates a chiral symmetry restoration with a temperature 
increase. In a narrow interval of temperatures 170 < T < 210 MeV the nucleon- 
u!-p and a excitation masses drop to zero (at T^a — 210 MeV) for the case 
ub = and similarly for finite n b ■ Masses of higher lying resonances also fall 
down but do not vanish at Tea- 

Thus at high temperatures and densities we deal with strongly interacting 
hadron matter which exhibits some (pre) critical properties. These properties 
can be related to the decrease of the a effective mass. The a excitation mass 
attains zero with the temperature increase till Tea, while it never decreases 
below a half of the bare mass at T = 0. The temperature dependence of 
the EoS at T near T^.^ reminds the "phase boundary" behavior, far away 
from the IG EoS. Particularities of the thermodynamic behavior near this 
"phase boundary" at high temperature were discussed in detail. In the case 
ub = the SHMC model results are compared with the lattice QCD data. 
For T > 170 MeV the pressure calculated in the SHMC model is much higher 
than that predicted by lattice simulations, that may give rise to problems in 
constructing a deconfinement transition if one wants to match this hadron EoS 
with a quark-gluon one at T > 170 MeV. The ways out of this problem have 
been discussed. The most natural way seems to introduce a different scaling 
law for baryon resonance coupling constants, as compared to those used for 
nucleons, resulting in an increase of the masses of baryon resonances in high 
temperature region compared to those used here. Then we could match our 
SHMC model EoS with that for quarks and gluons at quite high temperatures 
T > Tdec- In this case the quark liquid would masquerade as the hadron one. 
A similar idea has been recently discussed in application to hybrid stars in 
[83] . We demonstrated this idea by suppressing couplings of all baryons by 
factor 1/3, except nucleons. Then we easily match the lattice pressure up to 
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T ~ 240 MeV (Tea ~ 330 MeV for the case of suppressed couplings by factor 
1/3, except nucleons). 

The specific heat in SHMC model generates a sharp peak at T ~ 190 MeV 
simulating a strong crossover behavior. The peak is smeared out for suppressed 
couplings. It is interesting that at T > Tea the squared effective mass of the a 
excitation becomes negative. This can be interpreted as the occurrence of the 
instability with respect to a hot Bose condensation of a excitations. Similar 
possibility has been found in [21] in the framework of a different model and it 
was named as hot Bose condensation since the condensate, which stabilizes the 
system, appears for T > Tc rather than for T < Tc. Thus there arises a ques- 
tion on a competition between the hadron hot Bose condensation phase and 
the deconfined quark-gluon phase. All this may indicate that the description 
of the strongly interacting matter at high temperatures is more complicated 
than it was expected. In a more realistic model one should incorporate both 
hadron and quark-gluon degrees of freedom interacting with each other within 
a strongly correlated mixed quark-hadron state. Indeed, this idea is as old as 
twenty years [H] and recent quenched lattice QCD results show the existence 
of resonance structures above Tc. 

As an example for implementation of the SHMC model to HIC we examined 
the isentropic regime of this EoS which corresponds to expansion stage of 
nuclear system. A new method for extracting thermodynamic parameters of 
the freeze-out stage was considered using a shock-like model of freeze-out 
|77] . This method takes into account some elements of HIC dynamics and 
in-medium hadron modification. Obtained parameters prove to be a little bit 
different from those derived by the standard statistical model procedure for 
the energy range below the top AGS energy. The difference is within error 
bars for Tf^ and slightly above error bars for fifr extracted from data by 
the standard procedure. For more definite conclusions more precise data are 
needed as well as measurements of multistrange particles and inclusion of 
them in the scheme. Moreover, the results depend on the assumed model of 
the freeze-out. If we used the prompt freeze-out concept [22|76] we would 
obtain significantly larger difference between in-medium and ideal gas yields. 

If we want to apply the SHMC model for analyzing the particle ratios in the 
RHIC energy range, we should extend the used particle set to higher meson 
and baryon resonances. Here we don't do that since their coupling constants 
are unknown and we won't generate extra uncertainties. In addition, our aim 
is to demonstrate ability of the model in a broad (T, /i^) range rather than 
to study specific regimes. We are planning to return to this question in future 
publication. 

As we mentioned, if one moves to still higher temperatures, quark-gluon de- 
grees of freedom become deciding. We plan to match our hadron SHMC EoS 
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with quark-gluon EoS developing either a masquerade or a mixed phase sce- 
nario to cover a high T — fis range of the phase diagram. Then, we would like 
to insert this EoS into hydrodynamic codes to study dynamics of HIC. 

The SHMC model allows for arbitrary isotopic composition. It is interesting to 
apply it to collisions of isospin-asymmetric nuclei for studying isotopic effects 
and for describing initial stage of neutron star formation and cooling process 
(at temperatures T < 50 MeV and high baryon densities). 

The model allows one a further development. It is well known [T8|22] that pions 
and kaons have strong nNN, irNA, KNA, KNT,* interactions in the p-wave 
that we disregarded including only s-wave terms. These interactions can be 
incorporated into the scheme. For the sake of simplicity, the gas of quasipari- 
cle excitations was treated as a non-interacting gas. Feasibly this assumption 
is unrealistic when we deal with high temperatures. In the paper body we 
showed how one can generalize the model to the Hartree level. This can be 
done straightforward. However using RMF models and their generalizations 
we should always balance between a realistic and a practically tractable de- 
scriptions. Thus we postpone with further generalizations of the SHMC model 
to future work. 

Results are sensitive to the concept of freeze-out and it would be attractive to 
probe different freeze-out models. 

We used spatially homogeneous solutions. However initial equations of motion 
for mean fields are written in coordinate space. Thus in principle model al- 
lows to describe structures of possible mixed phases that may arise in neutron 
stars (cf. [29]) and feasibly in heavy ion collisions. With spatially inhomoge- 
neous solutions at hand one could apply the model to study atomic nuclei and 
confront it with many new experimental constraints. 

These interesting questions are however beyond the scope of the present paper 
and need further investigation. 
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Appendix A. Energy density of the gas of boson excitations 

Here we present expressions for the partial contributions to the energy density 
of the gas of boson excitations, Eq. (142|) . 

The energy density of the a meson excitations on the ground of the a mean 
field is 



oo 







Er'= UMp)). (64) 



with the dispersion relation 



c^<x(p) = V(mr"*)2 +p\ (65) 
and the effective mass (HSl): 



m% 7]^ dp 



(66) 



where Emf[/, o^o(/), i?o(/), ^] is given by Eqs. ([24D, ([25]), ([26]). All derivatives 
are taken here over variable / with fixed particle occupation. 

The energy density of the vector cu-meson excitations on the ground of the u 
mean field is 



J 2 

E^' = i2s. + 1) I ^ cuUp) f.Mp)) , (67) 



with the uj spin s^j = I and 

coUp) = v/(mr**)2 rnP-** = mj$^(/)| , (68) 

cf. Eq. (1441) . Similarly the energy density of the vector and iso- vector p-meson 
excitations on the ground of the p mean field is 

EP'^'-* = Ep+ + Epo + E,- = (2s, + 1) (69) 
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dp p 



with the p spin Sp = 1 and 



c^P± (p) = TV±9p x', Ro + ^{mr'y + P\ 
uAp) = V^(^^^rTT7, m^-^'* = m,|<|.,(/)| 

Conditions 9^a;^ = 0, S^p^ = are fulfilled. 
The energy density of the pion gas is 







dp p"^ 



27r2 

where for charged and neutral pions we have, respectively 



{p) = tV ± gl^ ujo ± g*p^ Ro + yjmf +p^ , 

uj^o (p) = ^{mlY ml = m^- gl^a 

with (yf*^ = due to absence of the u ^ 2n decay. 
The energy density of the kaon gas is 

E^j^'^ = Ek+ + Eko + Ek- + Eko 

J 2 

dp p 

[l^K+ip) fK+{^K+{p)) +^K"{P) fK"i^Koip))] 



27r2 


/dpp^ 
{uJK-ip) fK-iuJK-{p)) +UJko{p) fKo{^Ko{p))] 





where for charged kaons 



i^K± (p) = tV± glK ^0 ± g*pK Ro + + P^ , 

mK = mK- g*K(^ 
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and similarly for neutral kaons 



'^KyRoip) = ±9^K ^0 ± g*pK Ro + ym*^ +p'^ , (79) 
The r] contribution to the energy density, cf. [30], is given by 



^1 2 



where 



m*^ = I m: 



f2 " u I / 1 ' ^ P2 

b&{b} Jt^ / \ he{b} J-^ 

and the total baryon scalar density is Z^befb} < '^b >= Z]6g{6} 
The Bose distributions of excitations are 



nf-. 



exp[(y<^+7-/ir)/r]-i' 
^^l = /ii + gi(/ich + V) - gr - Qral^Ro, (83) 

i G {b.ex} = a, cu, p+, p°, p"; 7r+, 7r°, vr"; A'+, if", r/; 

Here pj = pstr for strange particles K and K* and pj = — pstr for their antipar- 
ticles. Qi is the boson electric charge in proton charge units, and we again did 
a gauge shift of the V variable, Q™'^ = +1 for particle, Q™'^ = —1 for antipar- 
ticle and QJ^^ = for the neutral particles (with all zero charges including 
strangeness). We take g^^.,, = 0, g^^K* = QpK* = 0, g^ji = for z = u; and gpi = 
for i = p. For simplicity m|^, = mx*, rn*, = m^i , m* = are assumed due 
to absence of corresponding experimental data. 

Note that the inclusion of the p-wave pion and kaon terms can be easily 
done. For that one needs to replace ci;^(p) and ujk{p) to more complicated 
expressions which can be found in Refs. [T8ll22f23j . However, we disregard 
such an interaction since its inclusion is beyond the scope of our RMF based 
scheme. 



Appendix B. Condensation of excitations 

There are two possible types of condensations of excitations. One type is the 
Bose-Einstein condensation which might occur for particles with a conserved 
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charge described by a finite value of tlie corresponding chemical potential. 
The number of condensed particles is determined by the value of the con- 
served charge (electric and/or strangeness in our case) even if one ignores 
self-interaction of bosons. Another type of condensation might occur for neu- 
tral particles. For them chemical potential is zero. The stability in this case 
is achieved by the interaction of the given boson species with other particle 
species and by the self-interaction. 

Bose-Einstein condensation. Let us consider charged bosons (a charge 
means either electric charge or strangeness). The density of the gas of the 
Bose excitations of the given species i is determined by the integral 



n. = g^j^m. (84) 



where gi is the degeneracy factor and fi is defined by Eq. (l82l) . For T > 
we have ^Ji*i{T) < m*{T) and there is no Bose-Einstein condensation. The 
critical point of the Bose-Einstein condensation is determined by the condition 
/i*(Tci) = m*(Tci). For T < Td the integral fl8^ diverges. Thus for T < Td 
the excitation number density Ui acquires a condensate term 



ni = Qi J fiip) + rioi . (85) 



Assume that the non-linear self-interaction term in the Lagrangian density 
is omitted. Then the relation noj = 2m*(T)|0^^p connects the condensate 
density and the mean field (jf^ (for simplicity we assumed the presence of a 
single condensate field). Then the energy density acquires extra term 



5E= noi{n,,T)m*. (86) 

The condensate does not contribute to the pressure. Contribution to the en- 
tropy of the Bose-Einstein condensate is also zero. Since —^*ni = —m*ni for 
T < Tcii the thermodynamic consistency relation holds. 

Generalization to the case, when the self-interaction term in the Lagrangian 
density is included, can be done as follows. The relation between //* and the 
mean field is found from the equation of motion for the mean field. As 
we have mentioned we should present the given boson field as = + 0'. 
In the gas approximation for excitations we should keep only quadratic terms 
(oc |0'P) in the Lagrangian density. Consider again a single (quasi) Goldstone 
excitation. By averaging £q*(0') over the equilibrium state we find 
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5U =< >- A|0"V/2 + 2A|0=f < |0f >, (87) 



< \<h'\^ >= riT 



where n^^ is the boson scalar density (sum of tadpole diagrams) for the given 
species. 

Minimization of the total energy results in equations of motion for the classical 
field and excitations. For classical field we obtain: 



101' = i^/mo, (89) 

where 

{u;, + ii*f-m*^-2\nl' = ^, (90) 

and LJc is the critical frequency given by + A** = ''^c- replace Uc + fi* 

by fi*(Tc) and put A = we recover the Bose-Einstein condensation result 
discussed above. 

The spectrum of particles now becomes 

(u; + /i*)2 = m*2+p2 + 2A|0^f . (91) 
Thus as compared to calculations done in absence of the condensate, we should 



replace u = \J m*^ + — fi* hj u = yjm*^ + 2A|0'^^P + ~ f^* in equations 
for particle excitation energies and add the terms A|0'^'|^/2 and — A|0'^'|^/2 
to the energy density and pressure, respectively. For charged pions the chiral 
symmetry consideration provides the relation A = ml/ f^. 

Condensation of neutral bosons. For neutral bosons, if strangeness and 
electric charge are zero, condensations may also appear. In this case fi* = 0. 
The critical point is determined by the condition m*{Tci) = 0. The stability 
is achieved only due to the non-zero meson-meson self-interaction with the 
positive coupling constant Aj for 0^ term. For t] and ip we use free dispersion 
relations, thereby they never condense in the framework of our consideration. 

In spite of the fact that condensate of the a-excitations appears for T > Tea 
rather than for T < Tea, it can be considered in the similar way as the cases 
discussed above. 
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